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We use the effective chiral Lagrangian to analyze the phase diagram of two-flavor twisted mass 
lattice QCD as a function of the normal and twisted masses, generalizing previous work for the 
untwisted theory. We first determine the chiral Lagrangian including discretization effects up to 
next-to- leading order (NLO) in a combined expansion in which / (Airf^) 2 ~ aA (a being the 
lattice spacing, and A = Aqcd)- We then focus on the region where m%/ '(4ir f n ) 2 ~ (aA) 2 , in which 
case competition between leading and NLO terms can lead to phase transitions. As for untwisted 
Wilson fermions, we find two possible phase diagrams, depending on the sign of a coefficient in the 
chiral Lagrangian. For one sign, there is an Aoki phase for pure Wilson fermions, with flavor and 
parity broken, but this is washed out into a crossover if the twisted mass is non-vanishing. For the 
other sign, there is a first order transition for pure Wilson fermions, and we find that this transition 

. extends into the twisted mass plane, ending with two symmetrical second order points at which 

the mass of the neutral pion vanishes. We provide graphs of the condensate and pion masses for 
both scenarios, and note a simple mathematical relation between them. These results may be of 

, importance to numerical simulations. 
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. I. INTRODUCTION AND CONCLUSION 

> 

There has been considerable recent interest in the twisted mass formulation of lattice QCD (tmLQCD). This for- 
mulation has several significant advantages compared to "untwisted" Wilson fermions: first, the twisted mass provides 
an infra-red cut-off for small eigenvalues of the lattice Dirac operator and thus avoids so-called exceptional configu- 
rations [HQ 1 second, for maximal twisting (i.e. a purely twisted mass term) physical quantities are automatically 
0(a) improved p[; and, finally, calculations of weak matrix elements are considerably simplified It may thus 

serve, for practical simulations, as an intermediate formulation between improved staggered fermions (which share 
the above advantages, but have the disadvantage of taking the square- or fourth-root of the fermion determinant) and 

'"" . chiral lattice fermions (which are more computationally expensive) . 

We consider here simplest version of tmLQCD which, by construction, is equivalent in the continuum limit to QCD 
with two degenerate flavors. At non-zero lattice spacing, however, the flavor group is explicitly broken from SU (2) 

i-C ■ down to U{\). This flavor breaking is analogous to the "taste symmetry" breaking of staggered fermions that plays 
, a major role in practical simulations. It is thus important to study the impact of flavor breaking in simulations 
of tmLQCD. This can be done analytically using the chiral effective theory including the effects of discretization. 
The methodology for doing so for Wilson fermions was worked out in Ref. 0. Here we provide the generalization 
to tmLQCD, determining the chiral effective Lagrangian to next-to-leading order (NLO). This extends the work of 
Ref. by including terms proportional to a 2 (a being the lattice spacing). 

There are many possible applications of the resulting effective Lagrangian. In this paper we focus on the phase 
structure. We extend the analysis of Ref. (and the similar considerations of Ref. 9J) into the twisted mass plane. 
For untwisted Wilson fermions, Ref. found two possibilities as the (untwisted) quark mass approaches the critical 
mass. In the first, the pion masses vanishes, and the theory enters an Aoki phase, in which flavor and parity are 
spontaneously broken. This is the scenario proposed long ago by Aoki 0, 0, , and supported by results from 
quenched simulations. The phase has width Am ~ a 2 , in terms of the physical quark mass m. The second possibility 
is that there is a first order transition at which the condensate flips sign, and the pion mass reaches a minimum, but 
is non- vanishing. In this scenario flavor and parity are unbroken. The choice of scenario is determined by the sign of a 
particular coefficient in the chiral Lagrangian [the coefficient is c-i defined in eqs. I|18ll9|l ]. the value of which depends 
on the gauge action and the coupling constant. Indeed, numerical work with unquenched Wilson fermions 01 finds 
evidence for an Aoki phase at strong coupling (consistent with the expectations of analytic results in the strong 
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(a) Phase diagram for C2 > (b) Phase diagram for C2 < 

FIG. 1: Phase diagram of tmLQCD. a and are proportional to the untwisted and twisted mass, respectively, in units 
proportional to the lattice spacing squared [see eqs. (I28l29l ll. The sign of the coefficient C2 determines whether flavor symmetry 
is spontaneously broken in the standard Wilson theory. The solid lines are first order transitions across which the condensate 
is discontinuous, with second-order endpoints. Figures 12141 show the dependence of the condensate and pion masses along the 
horizontal dashed lines. See text for further discussion. 

coupling, large N c limit 10]), but also finds that the phase disappears as the coupling is weakened. This suggests that 
the second scenario applies for moderate and weak couplings, at least with the Wilson gauge action. This conclusion 
is supported by the very recent work of Ref . |14| , who find clear evidence of a first order transition along the Wilson 
axis. 

The extension of the analysis of Ref. into the twisted mass plane turns out to be relatively straightforward. The 
twisted mass introduces only a single additional term in the potential, one that attempts to twist the condensate in 
the direction of the full mass. We find that any non-zero value for the twisted mass washes out the Aoki phase of the 
first scenario, as expected when the would-be spontaneously broken symmetry is explicitly broken. Thus the Aoki 
phase itself is confined to a short segment of the untwisted axis, as shown in Fig.^a). Traversing this segment in the 
twisted direction there is a first order transition, with a discontinuity in the chiral condensate. More interesting, and, 
in light of the results of simulations noted above, probably more relevant, is the impact of twisting on the scenario 
with a first order transition. Here we find that the transition extends a distance of size A/i ~ a 2 (where /i is the 
physical twisted mass) into the plane, as shown in Fig. ^b). The transition weakens as \/i\ increases, ending with 
second-order points where the neutral pion mass vanishes. 

Looking at the phase diagram, the two scenarios appear related by a 90° rotation. In fact, this relation is exact for 
the condensate and the neutral pion mass, aside from corrections of 0(a 3 ) in the chiral effective theory. In other words, 
approaching the critical mass along the twisted mass axis the neutral pion mass has exactly the same dependence 
on the twisted mass as it does on the normal mass in an untwisted Wilson simulation with the opposite value of c^. 
Similarly, the condensates in the two cases are simply related by a 90° twist. This identity does not, however, hold 
for the charged pion masses. These vanish in the Aoki phase, but do not vanish along the first order transition line 
for C2 < 0. Nevertheless, there is still a simple relation between the two scenarios: the charged pion mass-squareds for 
C2 < are obtained from those for C2 > by rotating by 90° and then adding a constant positive offset (proportional 
to \c 2 \a 2 ). 

The possible presence of phase structure extending in the twisted mass direction may hinder numerical simulations — 
there would be premature critical slowing down due to the second-order endpoint, and metastability beyond O nc 
can reduce these effects either by working at weaker coupling (since A/i oc a 2 ), or by attempting to find a gauge action 
which gives rise to a smaller value of C2 . 

The rest of this paper is organized as follows. In the next section, we present the two-step process of obtaining an 
effective chiral theory that describes the long distance physics of the underlying lattice theory of tmLQCD including 
discretization effects. We then, in Sec. lIIII explain how to use the resulting chiral Lagrangian to investigate the phase 
diagram of tmLQCD. 

An partial account of the work pres ented here was given in Ref. [l5l ] . Results on the masses and decay constants of 
pions will be presented elsewhere [T(j . 

Similar conclusions have been reached independently by Miinster 0, and Scorzato • 
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II. EFFECTIVE CHIRAL THEORY 



The theory we consider here is tmLQCD with a doublet of degenerate quarks. 2 In this section, we start by briefly 
reviewing the symmetry properties of the tmLQCD action. We obtain the effective chiral theory for pions by first 
determining the effective continuum Lagrangian at the quark level and then matching it onto the chiral Lagrangian. 



A. Twisted mass lattice QCD 

The fermionic part of the Euclidean lattice action has the form 0, 0] 

Sf = « 4 E ^(*) [I E 7/*(V£ + v M ) - r~ E K V » + m ° + ^nHo] rf>i(x), (1) 

where tpi and are the bare lattice fields (with "P standing for lattice and not indicating left-handed), and V p and 
V* are the usual covariant forward and backward lattice derivatives, respectively. The bare normal mass, m n , and the 
bare twisted mass, /Xo, are taken to be proportional to the identity matrix in flavor space. A possible "twist" in the 
Wilson term has not been included because it can be rotated away by appropriate changes in too and fig . This means 
that we are using the so-called "twisted basis" , which we find simplifies the subsequent analysis. We will not need to 
specify the irrelevant parameter r in our analysis, other than that it should be a non- vanishing constant, satisfying 
|r| < I for reflection positivity What we have in mind is the value r = 1 used in most numerical simulations. 

The action given in Q), combined with any hypercubically invariant gluon action, shares the following symmetries 
with standard Wilson fcrmions: invariance under gauge transformations, lattice rotations and translations, charge 
conjugation, and the U(l) transformations associated with fermion number. It differs from standard Wilson fermions 
in two important ways. First, the flavor SU(2) symmetry is broken explicitly by the fio term down to the U(l) 
subgroup with diagonal generator T3. Second, ordinary parity V is only a symmetry if combined either with a discrete 
flavor rotation 

U (x P ), U k (x)^ul(x P -ak), k = 1,2,3 
hon,2^i(xp) (2) 

-i-0;(a;p)Ti ! 27o 

[where xp = {— x, i), and our generators are normalized as rf — 1], or combined with a sign change of the twisted 
mass term 

V = V x {fio -» -Mo] ■ (3) 
The additional symmetries called 7^g P and 1Z$ x T>d in Ref. |^ will not be needed here. 




B. Effective continuum Lagrangian at the quark level 

Following the program of Symanzik pH |22| . the long distance properties of tmLQCD can be described by an 
effective continuum Lagrangian of the form: 

C eff = C Q + aC 1 +a 2 C 2 + --- ■ (4) 

The key constraint is that this Lagrangian need only be invariant under the symmetries of the lattice theory. We will 
work through the terms in turn, emphasizing the differences that are introduced by the twisted mass. The analysis 
parallels and extends that of Ref. Q • 

We consider first which consists of operators of dimension four or less, and is the Lagrangian which survives in 
the continuum limit. The lattice symmetries restrict its form to be 

£ = C g + 4>{Ip+ Z m — + «7 5 r 3 ^)V' - Z m —4>ii , (5) 
a a 



2 We do not consider the generalization to non-degenerate quarks discussed in Ref. |2(I . 
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where C g is the continuum gluon Lagrangian, and, as discussed below, (i oc /io- In particular, the parity-flavor 
symmetry Vp 2 requires parity and flavor to be violated in tandem: it forbids the flavor singlet parity violating terms 
tpjsip and F^Ffiu, as well as the flavor violating, parity even operator ipT^ip. The residual flavor U{\) symmetry 
forbids bilincars containing the flavor matrices r 12 - The coefficients Z m , \i and m c must be real in order to retain 
reflection positivity. Finally, the spurionic symmetry V implies that the the flavor-parity violating operator - 075 r 3' ! / ; 
comes with a coefficient odd in [io, here of linear order. 

The net effect is that the result J5J differs from that for Wilson fermions only by the /i term. In particular, the 
twisted mass is only multiplicatively and not additively renormalized Q . 

The relationship between the parameters in the continuum Lagrangian and those in the bare lattice Lagrangian 
has been given in Ref. [3(. The continuum fields ip are related to the bare lattice fields as follows: 

^-a- 3 / 2 Z[ 5 2 (a),ln(A rcg a)]^, (6) 

where Z is a matching factor, and A reg is the renormalization scale of the continuum theory. The physical quark mass 
is defined in the usual way, 

m = Z m (mo - m c )/a, (7) 
while the physical twisted mass is related to the bare parameters as 

/i = ^o/a=^pVo/a> (8) 

with Zp is the matching factor for the pseudoscalar density. Finally, the coupling constants in the continuum and 
lattice theories are related in the standard way. All the renormalization factors depend on g 2 (a) and, due to their 
anomalous dimensions, on ln(A reg a). We keep this (weak) dependence on the lattice spacing implicit. It should be 
borne in mind, however, that when, in the following, we refer to a power-law dependence on a, there will always be 
implicit sub- leading logarithmic corrections. We also note that the Z factors are mass independent (assuming that 
we are using a mass-independent regularization scheme in the continuum) and thus are the same for the twisted mass 
theory as for normal Wilson fermions. The same is true of rh c . 
In terms of the physical parameters, we finally have 

C = C g +iJj(]p+m + ij 5 T 3 n)ip . (9) 

In this theory, unlike on the lattice, the apparent flavor and parity breaking is fake, as fi can be rotated away by a 
non-anomalous axial rotation. This theory is thus equivalent to the usual, untwisted, two-flavor QCD with a mass 
m' = y/ m 2 + /i 2 . 

We next construct C\ . w hich contains dimension five operators. The enumeration of these operators is similar to 
that carried out in Ref. [23| , and we find 

Ci = biijjia^F^Tp + b 2 -4>(Ip+ m + i7sT 3 /i) 2 ^ 

+ b 3 rml}(p+ m + £75X3^)^ + b 4 mC g + b b m 2 '■ipip 

+ b 6 fnp{(p+ m + H5T3M), «75T3 ]"0 + brfx 2 ^ (10) 

The coefficients must be real for reflection positivity, and have a similar implicit weak dependence on a as do the 
Z factors discussed above. Note that the form of C\ is similar to that for untwisted Wilson fermions 0, except for 
the inclusion of the twisted mass in the operator (p+ m + «75T3/z), and the addition of the ^6 and 67 terms. 

A number of potential terms have been excluded by the V symmetry: m^iipip, m^ipi^^T^ip, /j^ipi'y^Tsip, ipD 2 ^^^ 
and , tpicr f _ lu F flv i^ 5 T3ip. The last of these, the "twisted Pauli term", requires a factor of fi and thus appears only in £ 2 - 
Charge conjugation symmetry enters at 0(a), forbidding a term similar to that with coefficient be but containing a 
commutator instead of an anticommutator. 

We can simplify C± by dropping terms that vanish by the leading order equation of motion, i.e. that which follows 
from Cq. This is equivalent to changing quark variables by an amount proportional to a. It removes the 62, 63 and 65 
terms. 3 



This is not strictly necessary. Terms vanishing by the equations of motion break the continuum symmetries in the same way as other 
terms which do not vanish. Since it is the symmetry breaking properties that matter when constructing the effective chiral Lagrangian, 
the form of the latter is the same whether or not the vanishing terms are kept. Nevertheless, we drop these terms as it simplifies the 
equations. 
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The situation with the 64, &5 and 07 terms is more subtle. On the one hand, they can be removed by 0(a) 
redefinitions of the parameters in £ (the coupling constant for 64, and the untwisted mass for 65 and 07). Note 
that for 67 this is an additive redefinition (m — > to + bjafi 2 ) which leads to a curvature in the critical mass m c as a 
function of the twisted mass. On the other hand, when one makes such a redefinition one loses direct contact with 
the underlying bare parameters, which are linearly related to <?, to and /i as explained above. If one wants to map out 
the phase diagram in terms of bare parameters, which is what one does in a numerical simulation, one should keep 
the 64, 65 and 67 terms. 

At this point, however, it is useful to anticipate the power-counting scheme that we use in our chiral effective theory. 
We consider to/A, ///A and aA as small quantities of the same size, and work to quadratic order in a combined chiral- 
continuum expansion. In this power-counting the three terms in question are all of cubic order, and can be dropped, 
avoiding the choice discussed in the previous paragraph. That the b 5 and 67 terms are of cubic order is manifest, since 
they involve two powers of quark masses in addition to the overall factor of a multiplying £\. The 64 term appears 
to be only of quadratic order, since it involves only a single power of the quark mass. However, this is a relative 
correction to a leading order term, and the leading order terms in the chiral effective theory arc themselves of linear 
order. Thus its contributions to the effective theory are of cubic order. 

Since we are working to quadratic order in the joint chiral-continuum expansion, we need to determine the form of 
£2, which contains dimension six operators. There are three such purely gluonic operators [24L l25j| . We do not write 
these down since, as explained in Ref. j2(J, they will lead to terms of too high order in our expansion. In particular, 
those which do not break the continuum rotation symmetry give rise to relative corrections proportional to a 2 . Since 
the leading term in the chiral-continuum expansion is of 0(m,a), these corrections are of cubic order and can be 
dropped. There are also gluonic operators which break the rotation symmetry down to its hypercubic subgroup, but 
these only enter the chiral Lagrangian at 0(a 2 m 2 ) because of the need to have four derivatives. 

The analysis of fermionic operators in the untwisted Wilson theory has been carried out in Refs. [53, ■ There 
are fermion bilinears and four-fermion operators. What is found, however, is that the only operators which break 
more symmetries than those already broken by £0 + £1 are the bilinears which break rotation symmetry down to its 
hypercubic subgroup. As for the gluonic operators discussed above, these contribute beyond quadratic order in our 
expansion. The contributions of the remaining operators to the effective chiral Lagrangian are of exactly the same 
form as those obtained by treating the effects at £\ at second order. Thus we do not list these operators. 

This leaves the new dimension six operators induced by the twisted mass. These, however, are all of cubic or higher 
order in our expansion. This is because they have at least one factor of /i in addition to the a 2 common to all operators 
in £2- An important example is the twisted Pauli term discussed above. 

Thus we conclude that, for the purpose of constructing the effective chiral Lagrangian to next-to- leading order 
(NLO) we need only keep the following terms at the quark level 

£ eff = £ g + to + ij 5 T 3 fi)i> + biaipia^F^ip , (11) 

with the relations to bare parameters given in eqs. J7J and (JSJ. The only difference from the result for the 
standard Wilson theory given in Ref. is the addition of a twisted mass term. 

C. Effective chiral Lagrangian 

Following Q , the next step is to write down a generalization of the continuum chiral Lagrangian that includes the 
effects of the Pauli term. As already noted, we use the following power counting scheme: 4 . 

1 > {m, p 2 , a} > {m 2 , mp 2 , p 4 , am, ap 2 , a 2 } > {m 3 , . . . } . (12) 

Here m is a generic mass parameter that can be either the renormalized normal mass to, or the renormalized twisted 
mass [i. There is no O(l) term in the chiral expansion for pions. The leading order (LO) terms are of linear order 
in this expansion, and the NLO terms of quadratic order. We work to NLO, which is sufficient for the study of the 
phase diagram. We stress that we can use the expansion in other regimes, e.g. n » s or m ~ a 2 , except that we then 
need to drop certain terms which become of too high order. 



4 Factors of A needed to make quantities dimensionless are implicit from now on unless otherwise specified 
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The chiral Lagrangian is built from the SU (2) matrix- valued E held, which contains the relevant low-energy degrees 
of freedom. It transforms under the chiral group SU(2)l x SU(2)r as 

E — * LEi?t, L € SU(2) L , Re SU{2) R . (13) 

The vacuum expectation value of Eo breaks the chiral symmetry down to an SU{2) subgroup. The fluctuations 
around Eo correspond to the pseudo-scalar mesons (pions): 

E(x) = E exp|i^7r a (a;)T a //| = E [cos(tt(x)//) + i-k[x) ■ t sm(n(x) / f)] , (14) 

^ a=l J 

where / is the decay constant (normalized so that f n = 93 MeV). The norm and the unit vector of the pion fields are 
given by ir = y/ir ■ n = 7r a 7r a and fr a — Tr a /^- 

The chiral Lagrangian can be obtained from the quark Lagrangian, eq. i|ll|) . by a standard spurion analysis. We 
must introduce a spurion matrix A for the Pauli term, as well as the usual spurion X f°r the mass terms. Both 
transform in the same way as the E field. At the end of the analysis the spurions are set to their respective constant 
values 

X — > 2B Q (m + ir 3 fj,) = rh + ir 3 jl and A — >2W a = a, (15) 

where B and Wq are unknown dimcnsionful parameters, and we have defined useful quantities m, ft and a. Note 
that the only change caused by the presence of the twisted mass is the appearance of the fi term in the constant value 
of X . Other than this, the construction of the chiral Lagrangian is identical to that for untwisted Wilson fermions. 

Because of this simplification, we can read off the form of the chiral Lagrangian for tmLQCD from Ref. j^, in 
which the Lagrangian for untwisted Wilson fermions was worked out to quadratic order in our expansion. The only 
extension we make is to include sources for currents and densities. The left- and right-handed currents are introduced 
in the standard way by using the covariant derivative -D^E = c^E — iZ M E + iEr^, and the associated field strengths, 
e.g. = dpl v — d v l^ + i[lfj„l v ], and enforcing invariance under local chiral transformations j^l^ij. Sources for 
scalar and pseudoscalar densities are similarly included if we write x — 2i?o( s + ip), with s and p hermitian matrix 
fields. Although the sources play no role in the discussion of the phase diagram, we have used them to obtain forms 
for various matrix elements of phenomenological interest ^BJ . These results will be described elsewhere 0] . 

Putting these ingredients together the resulting effective chiral Lagrangian is (in Euclidean space) 



C x = ^Tv(D^D^) - ^Tr( X f E + tf x ) - ^Tr(A^ + EU) 

- LiTr^EZ^Et) 2 - L 2 Tr( J D Al EA,£ t )Tr( J D Al EA,£ t ) 
+ (L 4 + L 5 /2)Tr(D M EtD M E)Tr( X t E + tf x ) 

+ L 5 {Tr [(D M EtD M E)( X t E + tf x )] - Tr(D^ D^)Tr( X ^ + ^ X )/2} 

- (L 6 + Ls/2) [Tr( X f E + 2 - L 8 {Tr[( X ^ + ^ X f] - [Tr( X f E + tf x )] * /%} 
-L 7 [Tr( X tE-Et x )] 2 

+ ti 12 T>(£ pi/ D M EU„Et + R^D^DJL) + L 13 Tr(L^ER^S) 

+ [W A + W 5 /2)Tt(D^D^)Tt(A^ + EU) - (W 6 + W 8 /2)Tr( X t E + E^Tr^E + E f i) 
+ W 10 Tv(D^D^ + DpYtDpA) - (W^ + W£/2) [Tt(A^ + SU)] 2 

+ contact terms . (16) 

Here the Lj's are the standard Gasser-Leutwyler low-energy constants of continuum chiral perturbation theory, and 
the Wi's and W/'s unknown low-energy constants associated with discretization errors. For the latter constants we use 
the notation of Ref. [23, till. In fact, these constants related to discretization are the same as those for the untwisted 
Wilson theory, assuming the same gauge action, since the twisting enters only through the parameters in the spurions. 

In writing Ijltjfl we have used various simplifications that result from the fact that the flavor group is SU(2). 
In particular, because A is, in the end, proportional to an element of SU(2) with a real proportionality constant, 
and because we do not use it as a source, we can take A+E + EU to be proportional to the identity matrix, and 
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Ti^A^E — Ti^A) to vanish. Thus the following possible terms vanish for SU(2): 

+ W 5 {Tt[(D^D^)(A^ + E+i)] - Tt(D^D^)Tt(A^ + Sti)/2} 

- W 8 {Tr[(x f £ + Stx)(i f S + EU)] - Tr(x f S + E^Tr^E + eU)/ 2 } 

- WrTrCx+E - E^Tr^E - E+i) 

- W^{Tr[(itE + Eti) 2 ] - [TrO^E + EU)] 2 ^} 

- W^Tr^E-Eti)] 2 . (17) 

Similarly, the operators multiplied by L5, L7 and L 8 alone in l|16|) do not contribute once we set x to 2£?o (?n + it fi) , 
although they are needed when using % as a source. They do not contribute in our subsequent study of the phase 
diagram. 

Finally, we comment on the W±o term, which is present only because we include both discretization errors and 
external sources. This term can, in fact, be removed using the equations of motion, but we have found that it 
provides a useful diagnostic in computations of matrix elements. 

Our result l|ltj|) agrees with, and extends, the work of Ref. [8j], in which the chiral Lagrangian for tmLQCD including 
effects linear in a was determined. Our generalizations are to include the term of 0(a 2 ), to make the simplifications 
due to using the group SU(2), and to include the sources for currents. In comparing our result to that in Ref. [j| it 
should be noted that the result in that work is expressed in the physical basis, while we use the twisted basis. The 
two results are related by a non-anomalous axial rotation. 

III. ANALYSIS OF PHASE DIAGRAM 

To investigate the phase diagram of tmLQCD, we are interested in the vacuum state of the effective continuum 
chiral theory. To the order that we are working, the potential energy is 

V x =- |Tr(E + Et) + H [Tr(E + £t)] 2 + |Tr[i(S - Et)r 3 ] 

+ |1 {Tr[ ? (E - Et)r 3 ] f + f§Tr[i(£ - Et)r 3 ]Tr(E + ^), (18) 

where the explicit forms of the coefficients, and their sizes in our power-counting scheme, are 
ci = / (m + d) ~ to + a, 

c 2 = -8 [(2L 6 + L 8 )m 2 + {2W 6 + W 8 ) a m + {2Wq + W 8 ) a 2 ] - m 2 + am + a 2 , 
C3 = / 2 A ~ M= 

c 4 = -8(2L 6 + L 8 )fL 2 - [i 2 , 

c 5 = 16 [(2£ 6 + L s ) mfi + {W 6 + W 8 /2) afi}~ fi{m + a) . (19) 

Note that there are no relations between the Ci coefficients: all five are independent at non-zero lattice spacing. The 
extra terms introduced by twisting are those with coefficients c 3 , C4 and C5. 

To determine the pattern of symmetry breaking as a function of the coefficients Ci, we repeat the analysis of 0. 
We distinguish three regions of quark masses, each successively smaller by a factor of a: 

(i) Physical quark masses: 1 ^> m ~ /i ^S> a. In this case, both discretization errors and terms of quadratic order 
in masses can be neglected, so that only the c\ ~ m and c 3 ~ fi terms need be kept. In this case the symmetry 
breaking is as in the continuum. In particular, the condensate Eo lies in the direction of the full mass term, so that: 

^ = (m + ir 3f i) = (to + iT 3 fi) 

\J TO 2 + /i 2 \J TO 2 + fl 2 

The pions are degenerate, with masses 



,uz = ^/m 2 + fi 2 [1 + 0(a, to, fj,)] . (21) 



,2 



Note that this result holds also if either m > o or /1 > o. 



8 



(ii) Significant discretization errors: 1 ^> to ~ fi ~ a. This is the parameter region for which our expansion is 
most natural. The LO coefficients Ci 3 still dominate over the NLO coefficients 02,45, but the 0(a) term in c\ cannot 
be ignored. Thus the results for region (i) still hold, except that the untwisted quark mass to must be replaced by 
the shifted mass ml — to + a, or cquivalently to — > to' = m + oWq/Bq. In other words, there is an 0(a) shift in the 
critical mass, as noted in Ref. 0. In terms of this new mass the coefficients become 

ci ~ to', C2 ~ to' 2 + am' + a 2 , C3 ~ //, C4 ~ /x 2 , C5 ~ fx(m' + a). (22) 

(iii) Aoki region: m' ~ ~ a 2 . In this region the nominally NLO coefficient C2 is of the same size as the LO 
coefficients €1,3; all are of 0(a 2 ). The other NLO coefficients remain suppressed by at least one power of a: C4 ~ a 4 , 
C5 ~ a 3 . The fact that LO and some NLO coefficients are comparable might suggest a breakdown in convergence. 
This is not so, however. The largest NNLO term is this parameter region is ~ a 3 , and thus suppressed by one power of 
a. Loop corrections are also suppressed, since they are quadratic in ml and \i (up to logarithms) and thus ~ a 4 . The 
key result that allows this reordering of the series is that the leading order discretization error has exactly the same 
form as the term proportional to to and so can be completely absorbed into to', to all orders in the chiral expansion. 

As in the untwisted Wilson theory, it is the competition between the LO terms and the NLO C2 term that can lead 
to interesting phase structure. Note that both ml and \i must be of 0{a 2 ) in order for such competition to occur; if 
either ml or fi is of 0(a) then one is in the continuum-like region (ii). 

To determine the condensate we must minimize the potential energy 1)18(1 . We parameterize the chiral field in the 
standard way: E = A + iB • r with real A and B satisfying A 2 + B 2 = 1, so that A, Bi € [—1,1]. Similarly, the 
condensate (the value of E that minimizes the potential) is written Eo = A m + iH m ■ t. In the Aoki region the 
potential is 5 

V x = -cxA - C3B3 + c 2 A 2 . (23) 

Although c\ = f 2 m' and C3 = f 2 fi can both take either sign, we need only consider the case when both are positive. 
The other possibilities can be obtained using the symmetries of V x : if c\ — * —c\ with 02,3 fixed (i.e. when ml — > —to'), 
the condensate changes as A rn — * — A m , i?3 jm — > i?3, m ; while if C3 — > — C3 with C1.2 fixed (i.e. when /i — * —fi), then 
A m — > A m , i?3, m — * —B3 >m . We do not need to specify the transformation of the other two components of the 
condensate since they vanish, as we now show. 
To do so it is useful to define r and 9 by 

A = rcos8, B 3 =rsm9, A 2 + B 2 = r 2 , B\ + B\ = \-t 2 , 0<r<l, (24) 

in terms of which the potential becomes 

V x = -r(ci cos 9 + C3 sin 9) + r 2 c 2 cos 2 6» . (25) 

We will first minimize this at fixed r, and then minimize with respect to r. Based on the symmetries just discussed, 
we consider only ci^ > 0. At fixed r, the term linear in r, — (ci cos 9 + C3 sin#), has its minimum in the quadrant 
< 9 < 7r/2, while the term quadratic in r, c 2 cos 2 9, has its minima at ±7r/2 for C2 > and 0, ir for C2 < 0. It follows 
that, irrespective of the sign of c%, the minimizing angle, 9 (r), also lies in the first quadrant: < 9 (r) < tt/2. Its 
actual value is given by the appropriate solution to 

CisinfloM - c 3 cos6'o(?') = c 2 r sm[26o(r)] , (26) 

and we denote the value of the potential at the minimum by V x>m i n {r). To minimize with respect to r we evaluate 
the derivative: 

dV x , m in(r) _ dV x 
dr dr 

where to obtain this simple form we have used eq. I)26[l. Since C3 > by assumption, and sin^o(^) > as 0o(r) lies 
in the first quadrant, we find that V xm i n (r) is a monotonically decreasing function of r. The symmetries discussed 



C3 



6=0n(r) 



(27) 



5 Henceforth, we will drop references to the 0(a 3 ) corrections. 
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above show that this holds for any c\ 3 . Thus the absolute minimum of V x is at r — 1, and the minimizing angle is 

e m = e Q (r = 1). 

At this point it is useful to introduce scaled variables 

c i / 2 ( m + o) 



|pa| 8\2W' + Wl\a 2 



m'/a 2 (28) 



In words, a is the shifted, untwisted mass in units proportional to a 2 , while (3 is the twisted mass in the same units. 
In the Aoki region a and (3 are of order unity. In terms of these variables, the equation to be solved to determine 
A m = cos Qyyi is 

a 2 + B 2 - 4 a 2 
Al T aA 3 m + 1 A 2 m ± aA m - — = 0. (30) 

[This is just eq. l|26|l for r = 1 after some manipulation.] Here the upper sign is for c 2 > 0, the lower for C2 < 0. One 
must pick solutions with A m real and satisfying — 1 < A m < 1, and of these choose that which minimizes V x . Given 
A m , the other components are given by 



Bx, m = B 2>m = , B 3>m = sign(c 3 ) v /r^4X, (31) 

where the sign of B 3m follows from the form of the potential, eq. (|23|1 . 

The pion masses are given by the quadratic fluctuations about the condensate, using eq. (|14|) . The unbroken U{1) 
flavor symmetry ensures that the charged pions are degenerate. We find 



dA m + c 3 B 3 _ m - 2c 2 A 2 m \c 2 \ (3 , 
m ni = m„ 2 = j2 = (32) 



Am„ = m W3 - m wi = — , (33) 



2c2B| 
f 2 

where to obtain the final form for charged pions we have used eq. J2HJ. It is straightforward to show that the none of 
the pion masses become negative for any values of the parameters. 

Before we show plots of the results, it is useful to discuss how the sign of c 2 affects the solutions. Recall that the 
sign of C2 depends on the gauge action, and is not known a priori. For the untwisted Wilson action, the sign has an 
important impact: the Aoki phase appears only for c 2 > 0, while there is a first order transition for c 2 < 0. Once we 
extend the theory into the full twisted mass plane, however, the two possibilities are related. We focus on the A — B 3 
plane (i.e. we set r — 1) since we know from above that the minimum of the potential lies in this plane for all c 2 . The 
potential in this plane is 

V x = —C\ cos 9 - c 3 sin 6 + c 2 cos 2 9 . (34) 

If we change variables as follows 

ci = c 3 , c 3 = -c[ , c 2 = -c' 2 , 6 = 6' - tt/2 , (35) 

so that cos 6 = sin 6' and sin 6 = — cosd' , then the potential has the same form 

V x = -c[ cos 6' - c' 3 sin 6' + c' 2 cos 2 6' - c' 2 , (36) 

aside from an overall shift. This implies that if we take the phase diagram for, say c 2 > 0, and rotate it anticlockwise by 
90°, we will obtain the phase diagram for c' 2 = —c 2 < 0, with the components of the condensate given by A' m = —B 3tm 
and £? 3 m = A m . The mass of tt 3 is given simply by its value at the rotated point, m'^ = m n3 . 

This discussion does not include fluctuations in the other two directions. To obtain the charged pion masses, 
however, we can use the general result eq. (|33|l for the splitting between charged and neutral pion masses, and the 
equality of the neutral masses, to obtain 

KJ^m^+Am 2 -Am'/ (37) 

2 , 2c 2 sin 2 6> m 2c' 2 sin 2 6' m 
— ''Sri t p p l-JOJ 

2 -2c' 2 cos 2 <4 2c' 2 sin 2 6' m , on , 
= m m + —2 ~ - p ( 39 ) 

= ml - 2 4 ■ (40) 
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A 




a 



FIG. 2: The global minimum of the potential, A m , as a function of a, for C2 > and (3 = 0, 1, 2, 3. 



mlf 2 /c 2 



a 







^^^^^ 6 


(3 = 3^^^ 










\ / 1 





a 



(a) Mass of ni and tT2 



(b) Mass of 7T3 



FIG. 3: Mass of the pions as a function of a, for C2 > and /3 = 0, 1, 2, 3. 



Thus we find a fixed offset (independent of c\ and C3) between the charged pion masses. The offset has the correct 
sign so that all pions have positive or zero mass-squared. In particular, if C2 > there is an Aoki phase with m ni = 0. 
This rotates (as shown in Fig. 1(b) I into a phase with — —2c' 2 /f 2 = +2c 2 / f 2 > 0. 



It is straightforward to solve the quartic equation and determine the condensate and pion masses. The only subtlety 
is that care must be taken when either [3 or a is set to zero. This can be seen, for example, from the final form for 
the charged pion masses, eq. <|32[) . which becomes 0/0 when (3 — ► 0. In either of these two limits, a — > or (3 — > 0, the 
quartic equation reduces to a quadratic. The solutions for (3 — > have been given in Ref. Q, and we do not repeat 
them here. Those for a — ► can be obtained from those for (3 — > using the 90° rotation in the mass plane just 
discussed. 

For the remainder of this article we illustrate the nature of the phase diagram by plotting the condensate and pion 
masses. We display results as a function of a for fixed (3 and C2. We only show results for (3 > since, as explained 
above, those for (3 < differ only by changing the sign of B 3m . In particular, the pion masses are symmetrical under 
(3 — > — f3. We could also consider only a > 0, but find it clearer to show results for the both signs. 

We begin with results for C2 > 0. Figure [21 shows the form of the identity component of the condensate, A m — 
Tr(Eo)/2, for = 0, 1, 2 and 3. The corresponding pion masses are shown in Fig. [31 In the untwisted theory ((3 = 0) 
there are second order transitions at a = ±2, as shown by the kinks in A m and the vanishing of the pion masses. 
The Aoki phase, with B^ m ^ 0, which breaks flavor and parity, and correspondingly has two Goldstone bosons, lies 
between these second-order points. Once [3 is non-vanishing, however, the transition is smoothed out into a crossover, 
and the pion masses are always non-zero. Flavor is broken for all a, with the charged pions heavier than the neutral 



vanishes and B 2 m 



— 1, so 



pion by 0(a 2 ), as given by eq. l|3lfl) . For the special case of a = (maximal twisting), A,- 
the mass-squared splitting is 2c2/ ' f 2 for all (3 (which becomes a difference of 2 in the units in the plots). 

It is now possible to pass through the Aoki phase by varying (3, i.e. by changing the sign of the twisted mass. When 
doing so there is a first-order phase transition, since B 3 ^ m jumps from — A^ n 



a/i -<*74 to -^T^AJ, 
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-V 1 -« 2 /4- 

We now consider the case where C2 is negative. Fig. 0] shows A m and the pion masses as a function of a for fixed 
values of (3. Fig. 4(a) and Fig. 4(b)| show the results for fi — 0. As discussed in Ref. 0, the condensate jumps from 



So = 1 (and thus A m = 1, B = 0) for a > to Eq = — 1 (and thus A m = —1, B = 0) for a < 0. This is a first order 
transition without flavor breaking, so all pions remain massive and degenerate. 

The remainder of Fig. 01 shows what happens at non-zero twisted mass. The effect of fi is to twist the condensate, 
so that there is a non-zero T3 component B^ m . There is, however, still a first order transition at which B^^ m flips sign 
between ±(1 — (3/2) (assuming (3 > 0). The neutral pion is now lighter than the charged pions due to the explicit 
flavor breaking. The neutral pion has a mass m 2 ^ = 2 1 C2 1 ( 1 — (3/2) 2 /f 2 at the transition, while, as noted above, the 
charged pions have a constant mass given by m 2 ri 2 = 2\c2\/ f 2 . The transition weakens as \(3\ increases, and ends 



with a second order transition point at (3 = ±2, at which the neutral pion is massless [see Fig. 4(f) |]. For larger (3 the 
transition is smoothed out. Note that, once away from the transition, for a = the mass-squared splitting between 
charged and neutral pions is 2|c2|// 2 . 

These plots illustrate the general result shown above, namely that the C2 < case can be obtained from that with 
C2 > by a 90° rotation and appropriate redefinitions. Indeed, the results for C2 > can alternatively be viewed as 
the plots for C2 < at fixed a with (3 varying, and vice versa, with the exception of the charged pion masses, which 
differ by a constant offset of 2|c 2 |// 2 . To illustrate this latter point we plot, in Fig. [SJ the pion masses for C2 < as 
a function of (3 for fixed values of a. Comparing to Fig. |3 we see the equality of the neutral pion masses and the 
constant offset in the charged pion masses. 
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(a) Global minimum, (i = 



(c) Global minimum, j3 = 1 




(e) Global minimum, j3 = 2 




(g) Global minimum, /3 = 3 
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(b) Pion masses, fi = 
m 2 / 2 /M 
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(d) Pion masses, /3 = 1 
m 2 / 2 /M 




(f) Pion masses, (i = 2 

m 2 / 2 /H 




(h) Pion masses, /3 = 3 



FIG. 4: Global minimum A m and pion masses as a function of a, for C2 < and /3 = 0, 1, 2, 3. The dashed lines are for 711,2 
and the solid lines are for 713 . 
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(a) Mass of 7ri and n2 (b) Mass of 7T3 

FIG. 5: Mass of the pions as a function of f3, for C2 < and a = 0, 1, 2, 3. 
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